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A theory of image-potential states is presented for the general case where these surface electronic 
states are resonant with a bulk continuum. The theory extends the multiple scattering approach of 
Echenique and Pendry into the strong coupling regime while retaining independence from specific 
forms of surface and bulk potentials. The theory predicts the existence of a well-resolved series of 
resonances for arbitrary coupling strengths. Surprisingly, distinct image-potential resonances are 
thus expected to exist on almost any metal surface, even in the limiting case of jellium. 
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1. Introduction 

At surfaces and interfaces of metals, atomic, molecular 
or other discrete electronic levels couple to a continuum 
of states in the volume. The accurate physical description 
of this coupling is central to understanding a wide variety 
of basic processes in catalysis, nanoscience or molecular 
electronics. Here, we consider one of the most simple, yet 
fundamental model systems of this kind, image-potential 
states, resonant with a structureless continuum. 

Experimental and theoretical studies of the ultrafast 
dynamics of electrons in image-potential states have 
vastly improved our understanding of electronic excita¬ 
tion and decay processes at surfaces of metals [T]. Elec¬ 
trons excited to these states experience the Coulombic 
image force perpendicular to the metal surface. For crys¬ 
tallographic faces which exhibit a gap of the projected 
bulk bands in the vicinity of the vacuum energy, this 
gives rise to a Rydberg series of states, characterized 
by hydrogen-like wavefunctions in the vacuum and ex¬ 
ponentially decaying Bloch waves in the bulk [2 : , ;3j. The 
conceptual simplicity of these electronic states and their 
well-defined properties have not only allowed to identify 
and quantify the important factors that govern their de¬ 
cay mu]. In fact, they have attained the role of a kind of 
drosophila of electron dynamics. As such, they serve as a 
benchmark system for the development of new ultrafast 
experimental techniques and as a reference for investiga¬ 
tions of more complex electron transfer processes [6HI5]. 

Recently, several experiments showed the existence of 
a well-defined series of similar states in the absence of a 
projected band gap QH EMU. Such image-potential res¬ 
onances are depopulated mainly by elastic electron trans¬ 
fer into the bulk whereas the classical image-potential 
states can only decay inelastically. Since resonant elas¬ 
tic channels are expected to dominate electron transfer 
processes at interfaces in most applications, model stud¬ 
ies of the decay of image-potential resonances are most 
interesting. These perspectives are, however, seriously 
hampered by the lack of a rigorous theoretical description 
which goes beyond appropriately tuned model potentials 
[151 [TiH . The multiple scattering theory of Echenique and 


Pendry |2j to free-electron-like metals as well as density 
functional theory at the level of the GW approximation 
predict just one broad resonance 23!21]. These results, 
which had been well accepted for many years, are now 
contradicted by recent experiments for Al(100) [151 . 

In this paper, we extend the theory of Ref. [2] in order 
to allow all scattering channels to interfere. It will be 
shown that such interference effects decisively change the 
solution in the regime of strong coupling leading to the 
full resolubility of the Rydberg series beyond the n = 1 
state. 


2. Theoretical Model 

We follow ideas first introduced by Feshbach for prob¬ 
lems in nuclear physics [22] and apply an open quantum 
system formalism. In this approach one considers the 
following non-hermitian effective Hamiltonian 

'H eS = n 0 -ivvK (i) 

Ho describes the unperturbed states. In practice, it is 
a N x N matrix with diagonal elements E n and all off- 
diagonal elements equal to zero (no configuration inter¬ 
action). V is & K x N matrix that describes the coupling 
to K continuum channels. The eigenvalues A of H Q s are 
complex, with E res = Re(A) denoting the maxima of the 
density of states and T res = — 2Im(A) their widths. 

This type of open quantum system formalism has suc¬ 
cessfully been used to solve problems in various fields of 
physics [231 [23]. Gauyacq and co-workers have applied 
it to interpret results they obtained for the interaction 
of two atomic helium levels as a function of the distance 
from an aluminum surface [25] and, more recently, to ex¬ 
plain the long excited-state lifetime of a metallic double 
chain adsorbed on Cu(lll) [26] . 

Before we apply the formalisms to the series of image- 
potential resonances, it is instructive to simply consider 
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two levels coupled to one continuum channel, 

«c=(“ f)- P) 

Without loss of generality the two levels are assumed to 
have energies = fE 0 /2, the lower one with stronger cou¬ 
pling and the upper one with a coupling that is smaller 
by a factor of /. The parameter a is a measure of the 
overall coupling strength of both levels to the continuum. 

One notices that the coupling matrix iVV t in Heg ([ 2 ]) 
does not only contain the diagonal elements ia and iaf 
describing a decay. It also contains the off-diagonal cou¬ 
pling terms ioty/J. As shown explicitly in the Appendix, 
the resulting interference is mediated by the continuum. 
In the coupled system, not only the two levels are affected 
by the continuum. Similarly, the initially structureless 
continuum is disturbed upon interaction with the two 
levels and this disturbance, in turn, acts back on the res¬ 
onances. Although l~L 0 s does not describe the continuum 
itself, it includes this interference effect in leading order. 

The complex eigenvalues of § are 

Ai,2 = T \yJ E l - oi 2 {l+f) 2 + i2aE 0 (l-f)-ia 1 -j^. (3) 

The dependence of this solution on the value of the 
coupling parameters a and / is illustrated in Fig. |T| 
One recognizes that the two resonances attract each 
other with increasing coupling strengths a. Whereas 
the width of the lower, more strongly coupled resonance 
Ti increases monotonously with a, the width T 2 of the 



FIG. 1. Energy positions Ei = Re(Ai) and half widths 
Fi/2 = — Im(Ai) of the two-level system (|2j [3| as function 
of the total coupling strength a(l + f) for Eq = 1 and / = 
|, |, 1. The dashed line denotes r tot /2 = a(l + /). In the 
special case f = 1 both resonances are centered at E = 0 for 

«(! + /) > I- 


more weakly coupled upper resonance attains a maxi¬ 
mum width around a(l + /) = 1. For this coupling, the 
combined half width of both resonances F tot /2 = a(l+f) 
equals the initial level spacing Eq. In the limit of large 
coupling, the lower resonance attains the total width 
a(l + /) whereas the other one becomes infinitely sharp 
(see appendix). This narrowing is known as resonance 
trapping [23 . 

The construction of the effective Hamiltonian |l]) for 
the series of image-potential resonances requires the en¬ 
ergy levels of the decoupled series and the coupling 
strengths to the metallic continuum. These values are 
obtained from the multiple scattering theory [2j. Within 
this approach a surface state is viewed as a wave trapped 
between the bulk crystal and the surface barrier, simi¬ 
lar to the modes of a Fabry-Perot interferometer. The 
repeated scattering at the crystal boundary and the sur¬ 
face barrier results in an amplitude 


1 - r B rc expi(4> B + 4>c) 

where </> B and 4>c denote the phase changes between in¬ 
cident and reflected waves at the surface barrier and 
at the crystal, respectively. r B and rc are the corre¬ 
sponding reflectivities. In the usual image-state problem 
r B = rc = 1 and the condition for bound states is 

<t>B + <fic = 27rn, n = 0,l,2,... (5) 


In the case of resonances, the loss of flux due to elastic 
electron transfer to bulk states results in rc < 1- This 
can be accounted for in terms of a complex phase (f>c = 
(f> c + i<f> c with an imaginary component </> c = — In rc [2] ■ 

Since this scattering model does not take a possible 
perturbation of the continuum by the resonances into 
account, a direct evaluation of eq. Q will not yield the 
proper energies and widths of the resonances when rc is 
small, i.e. in the case of strong coupling. In the weak cou¬ 
pling limit (rc —> 1), however, it can safely be employed 
to determine energies and relative coupling strengths, the 
only input needed to construct H 0 ff- In linear approxi¬ 
mation, the condition <®> is fulfilled for 

= <t>c ( 6 ) 

with r/2 being the energy dependent decay rate of the 
amplitude Q. 

For the energy dependence of the phase change <(> B at 
the image potential we use the WKB approximation m , 
<j)B — [(—8A 1 ) -1 / 2 — l]7r. Above the band gap 4> c = n (see 
below), and we obtain a Rydberg series of resonances 


E 


0 

n 


1 

32n 2 ’ 


n= 1 , 2 ,... 


(7) 


with a decay rate or full width at half maximum given 
by 


r u = 


-In rc(E°) 1 


n = 1, 2, 


27r 8 n 3 


( 8 ) 
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FIG. 2. Half widths and decay times of the first five reso¬ 
nances as function of coupling strength and reflectivity. Ar¬ 
rows indicate coupling strengths where maximum widths are 
reached. 



In the case of a simple metal it is sufficient to restrict 
oneself to one continuum ( I\ = 1) like in the above ex¬ 
ample of the two-level system. Since Eq. <[8j) must corre¬ 
spond to the solution of 0 in the limit of small coupling, 
the matrix elements V„ are determined by the condition 
\Vn \ 2 = r°/2. 

In addition to a n 3 -dependence, the T 3 depend weakly 
on the quantum number n via the energy dependence 
of the reflectivity rc- In order to facilitate a general 
discussion, we will neglect this weak dependence in the 
following. This approximation allows us to introduce the 
dimensionless parameter 

a=-(-lnr c ) (9) 

7r 


as a measure of the overall resonant coupling strength of 
the whole series and we obtain 


l^n | 2 = 


r° 

n 

2 


a 

32 n 3 ' 


( 10 ) 


Please note that in terms of the numerical matrix di- 
agonalization necessary to find the complex eigenvalues 
of eq. 0, the approximation of a constant rc, has no 
advantage and can easily be dropped in calculations for 
specific materials. However, it will become apparent be¬ 
low that for the most interesting cases of strong coupling 
the energy dependence of rc is indeed negligible over the 
limited range of energies E 3 . 


The resulting energies and widths of the first five res¬ 
onances are plotted in Figs. [2] and [3] as function of the 
coupling parameter a. Specific numerical values of inter¬ 
est are collected in Table [TJ In the limiting case rc —t 1 
ora->0 the eigenvalues of ([I) evidently give exactly the 
resonance spectrum 0 and 0. As expected from the 
above analysis of the two-level system, dramatic changes 
occur as soon as r 3 /2 + r 3 +1 /2 ~ T 3 = a/16 n 3 starts to 



FIG. 3. (a) False color contour plot of resonance spec¬ 

tra n = 1,.., 5 as function of coupling strength/reflectivity. 
For better visibility the Lorentzians of the individual reso¬ 
nances are superimposed with equal amplitudes (not areas); 
the dashed line indicates the maximum of n = 1 resonance, 
(b) Resonance spectra n = 2, ..,5 on an expanded scale; the 
transition between red and blue colors occurs approximately 
at half the maximum amplitude. 


approach the level spacing A 3 = E® +1 — E® ~ l/16n 3 , 
i.e. when a approaches unity j28) . 

In the regime of small coupling (a < 0.2) all resonances 
develop uniformly and simply broaden linearly as a func¬ 
tion of a (Fig. [3]). With increasing coupling the behavior 
of the n = 1 resonance and the other resonances become 
qualitatively different. Whereas the width of n = 1 in¬ 
creases monotonically and spreads over the whole spec¬ 
trum for a > 0.5, all other resonances reach a maximum 


TABLE I. Limiting values of energies A 3 , E 0° = E n (a—> 00 ), 
maximum widths r™ ax , minimum lifetimes r™ m and corre¬ 
sponding coupling parameter a(F™ ax ) for the first five reso¬ 
nances. 



E° n (eV) 

E™ (eV) 

rr x (eV) 

r-r (f S ) 

„ /pmax\ 

a(F„ ) 

n =1 

-0.8504 

-0.7363 

- 

- 

- 

n =2 

-0.2126 

-0.2967 

0.09545 

6.90 

0.6795 

n =3 

-0.0945 

-0.1221 

0.01920 

34.28 

0.5076 

n=4 

-0.0531 

-0.0593 

0.00662 

99.43 

0.4329 

71= 5 

-0.0340 

-0.0368 

0.00295 

223.10 

0.3871 













4 


width F™ ax and get narrower again for further increas¬ 
ing a. The higher the quantum number n, the sooner 
T™ ax is reached (Fig. [2]). The positions of the resonances 
also change as a function of coupling strength until they 
reach a limiting value Eff for large a. As illustrated 
by Fig. i the first resonance and the series of all other 
resonances effectively attract each other with increasing 
coupling. This attraction is larger, the lower the quan¬ 
tum number n. As a result, the resonances n = 2,3,... 
spread as a function of a. For the maximum coupling 
strength a = 2 plotted in Fig. [3j the energies E n have 
reached their limiting values already within 2% whereas 
the deviation from the uncoupled situation is more than 
30% for n = 2. 

Most importantly, if one compares the widths and the 
energy separation of the resonances n = 2,3,... (Fig. [3]), 
they are found to be well separated throughout the range 
of coupling strengths a = 0 —> oo. Even around a = 0.5 
where they reach their maximum widths T™ ax , the com¬ 
bined half-width of neighboring resonances stays well be¬ 
low their energy difference. This means - with the notable 
exception of the first resonance - that the integrity of the 
Rydberg series is not destroyed by strong coupling to a 
continuum. In contrast, stronger coupling even leads to 
a sharper spectrum of the higher resonances as compared 
to intermediate coupling. 

This, at first glance, surprising result is understood 
in a straightforward manner by considering the effective 
Hamiltonian ([T]) in the limit a —> oo. In that case the ma¬ 
trix VV^ determines the behavior of the system and one 
expects to obtain K states that take almost the whole 
coupling and N — K states that are almost decoupled. 
For our case of one continuum (Af=l), the most strongly 
coupled resonance (n=l) becomes the fast (open) chan¬ 
nel with Fi/2 ~ — Im{TV("Heff)} whereas n = 2,3,... 
become long lived (trapped) states. A reorganization of 
the system due to interference of the different decay chan¬ 
nel takes place already for intermediate values of a. As a 
consequence, the resonance spectrum never gets totally 
smeared out as one might intuitively expect. 

In terms of the basic physical processes, the trapping 
is nothing else than a back feeding of intensity from the 
distorted continuum into the more weakly coupled states. 
This is easily recognized in the example of the two level 
system discussed in the appendix. If the system in state 
|a) is coupled to the continuum at t = 0, the amplitude 
aft ) will start to decay. At the same time, the amplitude 
b(t), which was initially equal to zero, will start to in¬ 
crease in reaction to the distortion of the continuum (eq. 
Al). Since the amplitude bft) has the opposite phase as 
a(t), the coupling of state \b) to the continuum will in 
turn accelerate the decay of state |a). 


3. Results of Ag(lll), Al(lOO) and Al(lll) 

In our description, the only material-dependent pa¬ 
rameter is the reflectivity rc of a free electron wave at 


the surface. For surfaces of simple metals, like Ag(100), 
Al(100) and Al(lll), it can simply be obtained from the 
two-band model of the bulk electronic structure [25]. For 
a distance a of lattice planes along the surface normal 
and an sp-gap at the T-point of 2V g , the potential of the 
two-band model is 

V(z) = -Ho + V g e i9Z + V g e~ igz (11) 

with g = 2n/a, and Vo chosen such that the vacuum level 
is at zero energy. We restrict ourselves to a situation 
where the image-potential states are resonant with the 
upper band 

E{k) = -H 0 + (|) 2 + k 2 + yjv* + g 2 k 2 (12) 

The wave functions that solve the Schrodinger equation 
are Bloch waves in the crystal (z < 0) and plane waves 
in the vacuum (z > 0) 

'Fc = ae^- k)z + pe-W+Q*, z< 0 (13a) 

4> vac = Ae iqz + Be~ xqz , z > 0 (13b) 

with wave vectors 

k = \Je + Ho + (|) 2 — ^J~{E + Ho )<? 2 + V? 
q = \JE + Vq. 

The reflectivity rc is the ratio A/B of the in and out¬ 
going plane wave 'F vac in the vacuum, obtained from 
matching conditions at the interface (’F c /\FcU=o = 

'^vac/'^vac|z=o) [30J : 

_ A _ V g (q — fc + f) + V'(q — fc — §) 
r ° B V g (q + fc — f) + V'(q + fc + |) 1 ’ 

with 

V = Vg^ =gk+ \jv g + g 2 k 2 . 

Parameters used to describe the three surfaces 
Ag(lll), Al(100), and Al(lll) and the resulting reflec¬ 
tivities rc and coupling parameters a = —ir~ 1 lnrc for 



FIG. 4. Energy dependence of the reflectivity rc for Al(lll), 
Al(100), and Ag(lll) as calculated from the two-band model. 
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FIG. 5. Comparison of measured lifetimes (full circles) with calculated values (solid line) using reflectivities rc from the 
two-band model (Table [lit f° r Ag(lll) (a) and Al(100) (b); a common coupling parameter of a = 0.92 was used for Al(100), 
whereas the energy dependence of rc was taken into account for Ag(lll). For this surface, n = 1 (open circle) is an inelastically 
decaying gap state m- 


E = E vac = 0 are collected in Table |IT| Explicit values of 
the calculated resonance energies and lifetimes are given 
in Table |III| rc is equal to one at the edge of the upper 
band (k = 0) and decreases monotonously as the energy 
increases. This decrease is more rapid, the smaller the 
value of V g (Fig. [4]). 

Experimental data for the series of image-potential res¬ 
onances are available for Ag(lll) [T7j and Al(100) [T6| . 
Fig- E shows that the present theory is able to describe 
the measured lifetimes within a factor of two. Appar¬ 
ently, the two-band model predicts somewhat too small 
values of rc- In the case of Ag(lll), this results in 
shorter lifetimes than the experimental ones. In the case 
of Al(100), on the contrary, too small values of rc overes¬ 
timate the lifetimes of the high-n states that are already 
effectively trapped. 

For Ag(lll) n = 1 is still a gap state, rc is 50% near 
the vacuum energy. The corresponding coupling param¬ 
eter is a = 0.22. Trapping will thus be effective only for 
resonances with higher quantum numbers. For Al(100) 
the reflectivity is only 5% and the coupling parameter 
is a = 0.92. The theory predicts a very broad first 
resonance and widths of the high-n resonances that are 
substantially reduced as compared to the model of in¬ 
dependent decay (Fig. [2]). The quantitative comparison 


TABLE II. Reflectivities at the Vacuum Energy rc (E = 0) and 
coupling strengths a for Ag(lll), Al(100) and Al(lll); pa¬ 
rameters have been derived from the values given in Ref. m 
(band gap: E gap = 2V g , upper band edge: E up = E(k = 0)). 

Surface a (A) V g (eV) V 0 (eV) E up (eV) r c a 
Ag(lll) 4.085/\/3 2.150 9.3619 -0.4527 0.51 0.22 

Al(100) 4.049/^2 0.840 10.9723 -5.5450 0.054 0.92 

Al(lll) 4.049/^3 0.125 15.6485 -8.6425 0.0059 1.63 


with the experimental data of Refs. m and m shows 
good quantitative agreement for both surfaces, particu¬ 
larly if one considers that we neglect other decay chan¬ 
nels and make use of the most simple estimate for rc- 
Al(lll) is a surface very close to the idealized jellium 
model with a very low electron reflectivity (rc = 0.6%). 
With a resulting coupling parameter a = 1.63 the present 
theory predicts high-n resonances that are sharper and 
more long-lived than those of Al(100) and even those of 
Ag(lll). Electrons in these resonances are effectively 
trapped (Fig.[2]). 

It must be pointed out that the result for Al(lll) 
stands in marked contradiction with some of the es¬ 
tablished knowledge from the literature. Previous ex¬ 
perimental and theoretical work agree on the pres¬ 
ence of just one structureless feature near the vacuum 
level for Al(lll) arising from image-potential resonances 
gnuHEHSi]. However, the experimental resolution of 
Refs. m\m was less than 200 meV and not sufficient 
to resolve high-n resonances. A reinvestigation of well- 
defined Al(lll) with improved experimental possibilities 
would therefore be highly desirable. The Al(lll) surface 


TABLE III. Energies E n and decay times r n of the first seven 
image-potential resonances of Ag(lll), Al(100) and Al(lll). 



Ag(lll) 
E n (eV) r n (fs) 

E n (eV) 

Al(100) 

T n (fs) 

E n (eV) 

Al(lll) 
Tn (fs) 

n= 1 

- 

- 

-0.767 

0.4 

-0.745 

0.2 

n =2 

-0.211 

16.2 

-0.270 

7.3 

-0.288 

10.8 

n =3 

-0.095 

49.6 

-0.109 

41.6 

- 0.111 

65.9 

n=A 

-0.054 

116.2 

-0.058 

130.7 

-0.059 

221.4 

ii= 5 

-0.034 

229.2 

-0.037 

307.3 

-0.037 

502.9 

n =6 

-0.024 

402.6 

-0.025 

608.6 

-0.025 

1003.6 

n=7 

-0.018 

651.0 

-0.018 

1075.2 

-0.018 

1782.9 
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would also provide interesting opportunities to decouple 
the resonances by adlayers. The present theory clearly 
predicts that moderate decoupling should lead to shorter 
lifetimes of the high-n resonances, and not to longer ones 
as for image-potential states at gaps [551 . 


5. Conclusions 


In conclusion, we have shown that electrons in energet¬ 
ically adjacent states, strongly coupled to a continuum at 
a metal interface, do not delocalize independently from 
each other. Interference effects can lead to extremely 
long lifetimes of some states while others experience an 
accelerated decay. Criterium for strong coupling is a sig¬ 
nificant overlap of the broadened levels, not the absolute 
coupling strength. One expects it to be fulfilled in a va¬ 
riety of situations in nanoscience, catalysis or molecular 
electronics, e.g. when small clusters or organic molecules 
with many close lying levels interact with a substrate. 
As for the Rydberg series of image-potential resonances, 
our theory strengthens the prediction of Eclienique and 
Pendry [2], who stated that the whole series of states 
would be observable whenever the first two ones are sep¬ 
arated. With the notable exception of the first resonance, 
a series of clearly separated resonances should exist for 
all coupling strengths, i.e. on any well-defined metal sur¬ 
face. 
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We consider here only an interaction between the states 
|a) and \b) and the continuum states \/i) 

<■ a\V\n ) = V a , (b\V\n) = V b . 

This interaction is assumed to be energy independent; all 
other matrix elements (| Vj) vanish. The time dependence 
of the wave function of the full system 

I m) = a(t)\a) + b(t)\b) + £ c^e-^n) 

a* 

is determined by the Schrodinger equation = idt Ilk) 
which reads in matrix form 



/ a \ 


/E a V a ••• 14 \ 


/ a \ 


b 


E b V b ■■■ V b 


b 

id t 

c\e~ iuilt 

= 

V a 14 wi 


c\e~ iuit 


\c N e~ iu>Nt ) 


\V a 14 Wat/ 


\c N e~ iu>Nt J 


or as a set of coupled differential equations 


id t a(t) = E a a(t ) + 14 c M ( t)e lUfit (Ala) 

id t b(t) = E b b(t ) + V b ^2 c^e-^* (Alb) 

id t [c #1 (t)e _,u '' it l = V a a(t) +V b b{t) + w M c /i (t)e -,a, '‘ t . 

(Ale) 


We are only interested in the time dependence of the 
states | a) and | b) due to their interaction with the contin¬ 
uum, not in the evolution of the continuum states them¬ 
selves. We thus eliminate <4 from (Alai and (AlbI by 
integrating (AleI and obtain 


Cy.{t) 


dt'[V a a(t') + V b b(t')\e lul » t '. 


Appendix A: Effective Hamiltonian of the 
Two-Level System 


We further assume the continuum states \fi) to be equally 
spaced and to extend from —oo to +oo. Then the sums 
that appear in (|Ala|) and (lAlbh simply become 


In the literature, the effective Hamilton operator 'H e s 
for the general case of n levels interacting with several 
continua is usually obtained by means of the Feshbach 
projection operator formalism [2211221 [57] . Much of the 
essential physics, however, is already revealed by consid¬ 
ering two levels coupled to one structureless continuum 
[55] , In the following we give a simple derivation of "H e ff 
for the two-level system ([2]) and evaluate its properties. 
In this way, we hope to make the approach more easily 
accessible to the general surface science reader. 

The Hamilton operator of two electronic states |a) and 
| b) with energies E a and E b coupled to a continuum by 
an interaction V is 

H = Hq + V = E a \a)(a\ + E b \b)(b\ + |/i)(/t| + V. 

A* 


^c M (t)e- iw ^ = f dwc(w,f)e- iwt 

fl 

= -iJdu)J dt’[V a a(t') + V b b(t’)]e iui{t '- t) 
= -i2n[V a a(t) + 146(f)] 

and we arrive at two coupled differential equations for 
the amplitudes a(t) and b(t) 

(a(t)\ (E a - -i2nV a V b \ /a(f)\ 

* ‘ \b(t)) = ^ -i2nV a V b E b - i2nV b 2 ) ^6(t )) ' 

(A2) 

We identify the matrix in ( |A2[ ) with the effective Hamil¬ 
tonian for the two-level system |2]) with E a b = ^pE 0 /2, 
2nV 2 = a, and 2irV b = af. 
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One recognizes that the coupling to the continuum 
not only introduces a finite width of the levels, given 
by 27 tV _ 0 2 = a and 2 ttV 2 = af. In addition the two lev¬ 
els appear to interact with one another as described by 
the imaginary off-diagonal terms i2nV a Vb = ia\Jf. The 
derivation clearly reveals that this interaction is not a 
direct one, but one mediated by the continuum. The 
interaction with the levels |a) and | b) disturbs the con¬ 
tinuum states |/z). This time-dependent disturbance in 
turn acts back on |a) and |t>). This interaction becomes 
stronger with a stronger coupling to the continuum, and 
it is directly linked to the level widths. The derivation 
also shows that the effective Hamiltonian ([2]) will be a 
good approximation for the full Hamiltonian as far as 
the levels are concerned, when the continuum is struc¬ 
tureless and extends well beyond the energy spread of 
the broadened levels. This requirement is fulfilled in the 
case of the image-potential resonances considered in this 
work. 

The effective Hamiltonian ([2]) of the two-level system 
with energies ^fEo/2, overall coupling parameter a and 
relative coupling strength / (0 < / < 1) has two complex 
eigenvalues 


Ai, 2 = ~ « 2 (1+/) 2 +i2aE 0 (l-f) - 

([3} repeated) 

The resulting resonance energies and widths can be 
written in the form 

^ 1,2 = Re Ai _2 = T^pcos ^ (A3a) 

T 1j2 = -2ImAi j2 = a(l + /)±psin| (A3b) 

with p = p(E 0 , /, a) and ip = tp(E 0 , /, a). Obviously the 
total energy of the coupled system E to t = E\ + E 2 = 0 
does not change with the coupling parameter a while the 
total decay rate simply grows linearly with a 


coupling strength / 


A(a -S> 00 ) — E 0 + *^(^2 )■ 


For the degenerate case / = 1 the attraction is strictly 
quadratic until the two levels coincide for 2a > 1. While 
the system remains completely symmetrical in terms of 
the energies of the two resonances, the initial asymmetry 
of the half widths is seen to increase as a function of a. 

y ( a -> 0) = a + ^ 2/(1 — /) + ■•■ 

y (« -t 0) = fa - ^ fi 1 -/) + ■•• 


The width Ti of the more strongly coupled resonance 
increases monotonously. The width T 2 , however, reaches 
a maximum value ^max and goes to zero in the limit of 
infinite coupling. 


T lf x 

2 v ' 

r 2, ^ x 

— (a —> 00 ) 
2 v ; 


(1 + f)<* 


El f 

a (1 + /) 3 


El f 
a (1 + /) 


The maximum of T 2 is reached approximately when the 
real part of the argument under the square root in eq. (3) 
becomes zero. This is the case for 


a 


0 

max 


Eg 

1 + / 


(AS) 


i.e. for a value where the total decay rate T tot /2 = 
a max(l + /) is equal to the initial level spacing E 0 . For 
this coupling strength we have 


A («max) = Eo^j (jyy (A6a) 

r 2 « ax ) =Eo- A« ax ) (A6b) 


r to t = r 1 + r 2 = 2a(i + /). (A4) 

With increasing values of a, the energy separation of the 
two levels A = E 2 — E\ decreases. If one expands the 
solution into powers of a, one finds that the attraction is 
initially quadratic in a 


A (a 0) = E 0 + -^-a 2 - ... 

Eq 

and reaches a limiting value that depends on the relative 


These results concerning the maximum of F 2 are exact 
for / = 1. In this degenerate case, the widths of both 
levels increase strictly linear with a as long as 2a < 1. 
For 2a > 1, when their energies coincide, one of them 
increases faster, while the other one decreases and ap¬ 
proaches zero (s ee a lso Fig. 1 in the main text). For a 
ratio / < 1 eq. (A5) slightly overestimates the coupling 
parameter a for which the maximum of T 2 is reached. 
The error, however, is less than 12% and less than 4% in 
expression (A6b) which slightly underestimates r 2 !max . 
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